Emission Spectra of Fallback Disks Around Young Neutron Stars 
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ABSTRACT 

The nature of the energy source powering anomalous A-ray pulsars is uncertain. 



Proposed scenarios involve either an ultramagnetized neutron star, or accretion onto 
a neutron star. We consider the accretion model proposed recently by Chatterjee, 



Hernquist & Narayan, in which a disk is fed by fallback material following a supernova. 
We compute the optical, infrared, and submillimeter emission expected from such a 
disk, including both viscous dissipation and reradiation of A-ray flux impinging on the 
disk from the pulsar. We find that it is possible with current instruments to put serious 
constraints on this and on other accretion models of AXPs. Fallback disks could also 



be found around isolated radio pulsars and we compute the corresponding spectra. 



' We show that the excess emission in the R and I bands observed for the pulsar PSR 

0656+14 is broadly consistent with emission from a disk. 
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1. Introduction 



In the past few years, the realm of astronomical objects has been enriched with the discovery 
of about half a dozen of the so-called anomalous A-ray pulsars (AXPs) (Mereghetti & Stella 
1995; van Paradijs, Taam & van den Heuvel 1995). AXPs are sources of pulsed A-ray emission, 
with persistent luminosities L x ~ 10 35 — 10 36 erg s _1 and soft spectra. Their periods lie in a very 
narrow range, between 6 and 12 seconds and their characteristic ages are of order 10 3 — 10 5 years. 
To date, about half the known AXPs have been associated with Galactic supernova remnants. 
(For a recent review of the properties of AXPs, see Mereghetti 1999.) 

While AXPs are almost certainly neutron stars, their properties establish them as a distinct 
class from binary A-ray pulsars and radio pulsars. Compared to binary A-ray pulsars, AXPs 
have lower luminosities and exhibit a narrow distribution of periods. Unlike young radio pulsars, 
AXPs have rather longer periods and appear to be radio quiet. One key to understanding these 
differences is to try and identify the energy source which powers the A-ray emission. It is quite 
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clear that this energy cannot be provided by rotation (as it is in radio pulsars). For values of P and 
P that are characteristic of AXPs, the rate of loss of rotational energy is \E\ = 4tt 2 IP/P 3 « 10 32,5 
erg s _1 , orders of magnitude smaller than the observed X-ray luminosities. 

Two broad classes of models have emerged to account for the X-ray emission from AXPs. In 
one, AXPs are hypothesized to be isolated, ultramagnetized neutron stars, with field strengths 
in the range 10 14 — 10 15 G, a.k.a. "magnetars" (Duncan & Thompson 1992). Field strengths 
of this order are consistent with those inferred for soft gamma repeaters, on the basis of timing 
measurements, assuming that these objects lose rotational energy in a manner similar to radio 
pulsars (e.g. Kouveliotou et al. 1998, 1999). If the magnetar interpretation of AXPs is correct, 
the X-ray luminosities could be powered by either magnetic field decay (Thompson &: Duncan 
1996) or by residual thermal energy (Heyl & Hernquist 1997a, b). The observed X-ray luminosities 
place constraints on both mechanisms. If the emission is powered by residual thermal energy, then 
the envelope of the star must consist primarily of light elements such as hydrogen and helium 
(Heyl & Hernquist 1997b). On the other hand, if the emission is powered by magnetic field decay, 
then a value of B ^ 10 16 G is required, unless non standard decay processes are invoked (Heyl & 
Kulkarni 1998). 

The second class of models proposes that the X-ray emission from AXPs is powered by 
accretion, and unusually high values of the magnetic field are not required. Accretion can occur 
from binary companions of very low mass (Mereghetti & Stella 1995), from the debris of a 
disrupted high mass companion (van Paradijs, Taam & van den Heuvel 1995; Gosh, Angelini 
& White 1997) or from the interstellar medium (Wang 1997). Each of these three possibilities 
has been challenged by observational and theoretical considerations. Searches to detect binary 
companions of AXPs have failed thus far, placing severe constraints on companion masses (e.g. 
Wilson et al. 1998; Mereghetti, Israel & Stella 1998). If the emission is powered by accretion 
from a disrupted binary companion, it is not clear why AXPs should be associated with young 
supernova remnants. Finally, for the specific object (RXJ 0720.4-3125) considered by Wang 
(1997), Heyl & Hernquist (1998) have demonstrated that the X-ray luminosity can be explained 
by cooling, without requiring that this neutron star be old and have experienced magnetic field 
decay. 

Recently, Chatterjee, Hernquist & Narayan (1999; hereafter, CHN) have studied another 
accretion model, in which neutron stars accrete from disks which formed from fallback material 
after a supernova explosion. (Similar proposals have been studied recently by Alpar [1999] and 
Marsden et al. [1999].) The possibility that material might fall back onto a neutron star following 
a supernova explosion and settle in a disk is not new (e.g. Woosley 1988, Chevalier 1989). In 
particular, Lin, Woosley & Bodenheimer (1991) suggested that this process might account for the 
presence of planets around some radio pulsars. Later, in their discussion of accretion from the 
debris of a disrupted companion, van Paradijs et al. (1995) noted that fallback could provide 
material to form disks around young AXPs, but they did not analyze this possibility in detail. 
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In principle, all neutron stars could acquire debris through fallback, but subsequent 
accretion can occur only under specific circumstances, depending on the relative locations of the 
magnetospheric radius, the light cylinder radius, and the corotation radius. CHN related these 
conditions to physical parameters such as the initial mass of the disk, M<j, the initial period of 
the neutron star, Pq, and the strength of the magnetic field, B. CHN found that neutron stars 
would evolve into accreting sources only for sufficiently strong magnetic fields, while objects with 
weaker fields would be unaffected and would, presumably, become ordinary radio pulsars. For 
the particular values Pq = 0.015 seconds and = O.OO6M , for example, CHN found that the 
outcome depended on whether or not the field is greater or smaller than B = 4x 10 12 G. Namely, 
neutron stars with B > 4 x 10 12 G became AXPs, while those with B < 4 x 10 12 G became radio 
pulsars. Other choices for Pq and yield similar behavior, and lead to different values for the 
critical magnetic field separating AXPs and radio pulsars. 

In their analysis, CHN demonstrated that the time varying nature of accretion from a fallback 
disk provides a natural explanation for the fact that AXP periods increase with time. Indeed, if 
the accretion rate declines as M oc t~ 7 / 6 , close to the similarity solution derived by Cannizzo, 
Lee & Goodman (1990), then the rate of increase of the period will be identical to spindown by 
magnetic dipole radiation. Thus, timing analyses of steady spindown alone cannot distinguish 
between these two interpretations. Moreover, in contrast with the magnetar model, the time 
varying accretion studied by CHN relies on a standard population of isolated neutron stars having 
magnetic fields similar to those inferred for radio pulsars, and therefore it does not require the 
existence of a separate class of neutron stars. 

How can we distinguish between the time-dependent accretion interpretation and the 
magnetar model, given that steady spin-down is accounted for by both? One possibility may be 
by searching directly for disk emission, at wavelengths much longer than A-rays. Accretion disks 
dissipate energy by viscous processes at all radii, and some of this energy will be radiated at long 
wavelengths, especially from large radii. Moreover, the central A-ray emitting source irradiates the 
outer disk, and most of this energy is reradiated at longer wavelengths. In this paper, we compute 
the emission spectrum expected for accretion disks like those required to explain AXPs, and we 
apply the model of CHN to make predictions for the flux in various bands. The predicted flux is 
within the range of sensitivity of current instruments, and can be used, therefore, to directly test 
the accretion model. We note that while our specific estimates are made using the disk models of 
CHN, our arguments should apply in general to all accretion models for AXPs regardless of the 
source of the disk mass. Thus, sensitive observations at long wavelengths may provide powerful 
tests of all variants of the accretion hypothesis for AXPs. 

A prediction of the CHN analysis is that fallback disks could also be found around radio 
pulsars, although the properties of these disks are not well-constrained. Consequently, we apply 
our methods to calculate long wavelength emission from disks around radio pulsars. In particular, 
we consider the middle-aged radio pulsar PSR0656+14, which has emission in the R and I bands 
in excess of that expected by extrapolating the thermal blackbody spectrum observed at A-ray 
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wavelengths. We show that the excess emission is broadly consistent with that produced by an 
accretion disk. 



2. Model 

The evolution of the fallback material settling into a disk can be treated in a similar manner 
to the evolution of the accretion disk formed by material captured by a black hole after the tidal 
disruption of a star. Numerical simulations by Cannizzo, Lee & Goodman (1990; hereafter CLG) 
show that after an initial transient phase, the evolution of the disk becomes self-similar, with the 
accretion rate M decreasing as a power law. Exact self-similar solutions for the evolution of the 
disk were derived by Pringle (1974) assuming that v oc r p S 9 , where v is the viscosity and S the 
surface density. If the so-called a disk model is adopted (see CLG for details), then a solution is 
obtained when p = 1 and q = 2/3, yielding 

v = CrY?l' i , (1) 

with 

\[im v ) \l2acGMj 1 ' 

Here K es is the electron scattering opacity 0, M the mass of the neutron star, and a is the standard 
dimensionless viscosity constant that is expected to be ^ 1 on physical grounds (Shakura h 
Sunyaev 1973). We adopt the value a = 0.1 in our numerical calculations. The corresponding 
surface density evolves according to 



+ \ -15/26 

S(r,t) = S (-J / 



(3) 



where 

/(«) = (28)- 3 / 2 n- 3 / 5 (l - u^f' 2 . (4) 



The constants So, an d to are constrained by the condition 

,2/3 

The total mass of the disk decreases with time as 



M d (t) = (28)- 3 / 2 ^r6% (rV 3/16 • (6) 



7 u "Vto 



2 CLG considered more general opacity laws, such as the Kramers opacity, with n(p,T) — KopT~ 7 ^ 2 , but they 
found their results to be relatively insensitive to the particular expression they chose. Therefore, for simplicity we 
adopt a constant opacity, k cs . 
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The accretion rate is then given by M&it) = (—3/16) M&(t)/t. We define the initial mass of the 
disk to be = M^T^), where is of order the dynamical time of the disk early on. Following 
CHN, we select = 1 msec. Having chosen a value for rg, the constant So is then derived from 
equation (||) at t = T^. Using the relation (||), we then find 

So « 155 (M d ) 8 / 7 (T d Cf/ u r - 5/2 . (7) 

Notice that, with the choice (0), (i.e. So oc r 5 ^ 2 ), the self-similar solution for S, as given by 
Eq. (||), becomes independent of the numerical value of the parameter ro, and therefore is solely 
determined by M® and Td- However, as Equation (0) shows, the dependence on T d is much weaker 
than the dependence on M® 0. Following CHN, we keep the value of T d fixed and use only M® as 
a parameter. 

The solution (||) extends from u = to u = 1. However, an accretion disk can exist only 
for radii beyond the magnetospheric radius R rn ~ 6.6 x 1Q 1 m~ 2 /' 1 (Frank, King & Raine 
1992), where B\2 = B/(10 12 G) and rh = M/Me is the mass accretion rate in Eddington units. 
We therefore take the minimum radius of the disk to be r m j n = R m . Similarly, whereas the 
mathematical solution sets the disk outer boundary at u = 1 (corresponding to S = 0), in practice 
a disk can extend only up to a radius r max (S m i n ), where densities S <^ S m i n will be destroyed by the 
ram pressure from the interstellar medium. We assume that a disk can exist as long as the incident 
flux is not sufficient to significantly perturb the angular momentum of the matter orbiting in it. 
We find that this condition is satisfied for densities S ^ S m ; n «5x lCP 7 r 2 4 M[~J no v|oo g cm -2 
where ri4 = r/(10 14 cm), M1.4 = M/IAMq, no = nf (1 gm/cm 3 ) is the particle density in the 
interstellar medium, and V300 = f/300 km/sec. Therefore interstellar ram pressure affects the 
outer boundary of the disk only very marginally. 

At any position and time, the effective temperature of the disk due to viscous dissipation is 
given by (CLG) 



^ 9^(r,i)S(r,<)GM ^ 4 / 3M d (t)GM \ 
Tes(r,t)=[ ^3 ) o __ 3 ) 



Besides emission due to internal energy dissipation, a considerable amount of flux from the disk 
can arise as a result of reradiation of X-rays impinging from the central object. Vrtilek et al. 
(1990) derive an analytical expression for the effective temperature of an irradiated disk, under the 
assumption that the disk height h oc r n and irradiation is the dominant form of heating. If Lx is 
the X-ray luminosity of the central source, and u> = \J GM/r 3 is the Keplerian rotation frequency 
of the disk at position r, then n = 9/7, and the temperature profile is given by 



T x (r) 



1 UttctGM 



- 7 'JLY 17 f^i] 2/7 /R -^ 3/7 



^ 23200 (-4) (*®\ . (9) 



3 Large variations in the choice of Td would lead to the same numerical results with small adjustments in Mj. 
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The last equation is for M = IAMq. Here, / is a factor which contains the uncertainty in the 
disk structure (including the absorbed fraction of radiation impinging on the surface) . We adopt 
the estimate given by Vrtilek et al. (1990) of / ~ 0.5. In the cases that we consider, irradiation 
typically dominates the spectrum for Lx <; 10 34 erg s _1 ; for lower values of Lx (but generally 
higher than ~ 10 31 erg s _1 ), the inner parts of the disk (where the optical emission is mostly 
generated) are dominated by viscous heating, while irradiation dominates in the outer regions of 
the disk (where most of the longer- wavelength emission comes from). Therefore it is reasonable to 
keep the approximation (||) also in these cases and compute the contribution to the flux of the two 
terms separately. 

The total flux is obtained by integrating the emissivity over the entire surface of the disk: 

Y rrmax r2ir 

F = / drr d(j)[B^(T c s) + B u (Tx)] x [cos ^sin/3(r) sin a + cos (3{r) cos a] . (10) 

d J rmin J 

Here a is the angle between the line of sight and the normal to the disk at r = 0, and 
(3{r) = dh{r)/dr is the tilt angle of the surface of the disk, which we determine at every position 
in the disk from the equation of hydrostatic equilibrium in the z-direction. 

The opacity in disks made of fallback material is likely to be higher than in disks with cosmic 
abundance, due to the expected high metal content. Following Lin, Woosley & Bodenheimer 
(1991), we adopt an opacity of ~ 10 3 cm 2 g _1 in our calculations. Notice, however, that in our 
model the predicted spectrum is insensitive to the precise value of the opacity, as the temperature 
is independent of it. The size of the disk is slightly larger for high n es , but this has little influence 
on the final result. What is important, however, is that a high opacity ensures that the disk is 
optically thick at all radii, and therefore the emission can be computed as blackbody. A more 
significant uncertainty derives from the possible existence of thermally unstable regions in the 
disk. For disks with cosmic abundance, the unstable temperature ranges have been studied in 
detail (i.e. Cannizzo 1993 and references therein). For high metallicity disks the situation is likely 
to be quite different, as cooling depends on the abundances. 

3. Emission from disks around AXPs and radio pulsars 

In the model that we consider (CHN), AXPs and radio pulsars are assumed to be drawn 
from the same underlying population of neutron stars. In one example considered by CHN, for 
a fallback disk with mass MQ = 0.006Mq, and an initial period Po = 15 msec, the evolution of 
the accreting disk will lead to an AXP phase if the magnetic field of the neutron star is Byi ^ 4. 
Smaller values of the magnetic field would lead to the more common radio pulsar phase. In the 
following, we consider these two situations separately, and compute the emission spectrum from 
the disk at various ages. 
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3.1. AXPs 

Here we consider a disk of mass M° = 0.005M© at an angle of 60° with respect to the line of 
sight, and a magnetic field Byi = 8 for the neutron star. According to CHN, such objects will go 
through a phase of accretion-induced X-ray emission and will appear as AXPs. We describe the 
X-ray luminosity as 

m Le if m > 0.01 
100 m 2 if m< 0.01 



The rapid decline of the luminosity for m < 0.01 is to allow for the fact that the system is likely 
to switch to an ADAF phase, at which point the accretion luminosity Lx decreases faster than M 
(Narayan & Yi 1995; Esin et al 1997). Some studies have suggested, in fact, that in this phase 
much of the mass could be ejected prior to reaching the surface of the star (e.g. Blandford & 
Begelman 1999, Menou et al. 1999, Quataert & Narayan 1999). The transition radius between 
the ADAF and the disk is not well determined. Moreover, it is likely to depend on the metallicity 
of the gas in the disk, which is very far from the standard value. For simplicity, we keep the 
transition radius equal to the Alfven radius (which at later times is on the order of a few thousand 
Schwarzschild radii). We will discuss the consequences of this uncertainty on our results. 

Figure 1 shows the unextincted flux at different ages from an AXP at a distance of 5 kpc. 
The peak in the submillimeter range is due to the effect of irradiation. At early times (T <^ 10 4 
yr), irradiation dominates over viscous dissipation in the entire disk. At later times, when the 
X-ray flux is much lower, viscous heating becomes dominant in the inner parts of the disk, and 
this produces an enhancement in the emission in the optical region of the spectrum. The model 
predicts a flux of a few mJy at frequencies around 10 12 Hz, which is within the range of the 
SCUBA instrument. The infrared emission is also generally above the sensitivity limit of NICMOS 
on the Hubble Space Telescope. In the J band (1.1 /im), the limiting magnitude for a 5cr detection 
in 60 minutes of exposure is mag ~ 25, which corresponds to a flux of ~ 10~ 4 mJy = 10 -30 erg 
cm -2 s _1 Hz -1 . 

Presently, about half a dozen AXPs are known. In Table 1 we summarize the properties of the 
nearest ones, and compute the expected emission in various bands, assuming that the accretion 
model is correct. As an illustration, we take = O.OO5M0, B12 = 8, and an inclination angle of 
the disk i = 60° in all cases. X-ray luminosity, age and distance are taken from observations^], and 
extinction by dust has been taken into account based on the measured column density of hydrogen 
and an average extinction for the Galaxy^. Following Pei (1992), we fit the extinction curve by 

AW = k(W (ipIfO • (12) 



4 The observed X-ray luminosity is consistent with the value predicted by the accretion model of CHN. 

5 Needless to say, given the uncertainties in both theoretical and measured parameters, the numbers in the table 
must be simply considered as rough estimates 
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where 

6 



e(A)_ £(A/*)* + cW'+fc' (13) 



i- 



Here a,i,bi,rii, Aj are fitting parameters (that can be found in Pei 1992), and k is the dust-to-gas 
ratio. For the Galaxy, k = 0.78. 

Observations in some bands have already been made for a couple of the objects in the table. 
In the case of IE 2259, our predicted flux is below the upper limits set by Coe & Pighling (1998) 
in the J band (^ 19.6), and is somewhat higher than the limit of ~ 18 mag set in the K band. For 
the case of 4U 01412+615 reported by Mereghetti (1999), we get a flux comparable to the flux 
limit of 25 mag in the V band, and somewhat higher than the limit of 17 mag in the K band. We 
need to point out, however, that in our model the flux in the higher energy part of the spectrum 
is rather sensitive to the choice of r m ; n , which we take equal to the Alfven radius. If the inner 
bound of the disk is at a radius larger than this, then the optical and near infrared emission will 
be lower. The submillimeter flux, on the other hand, is produced at large radii, and therefore is 
quite insensitive to the specific value of r m - m . All fluxes are sensitive to M9, for which we have 
assumed the value 0.005 M suggested by CHN. It is unclear how much freedom is available in 
this parameter, but it is probably not very large. Another uncertainty is the inclination angle of 
the disk, for which we have assumed the value i = 60°. Furthermore, the distances to most of 
these objects are not well constrained. Given all these uncertainties, more and deeper observations 
in several bands are needed before drawing definitive conclusions. 

Observations in the optical (R and B band) have been made for an object not in the table, 
RXJ 0720.4-3125, which shares some characteristics with the known AXPs (Kulkarni & van 
Kerkwijk 1998), such as a long period (8.39 sec), but which has a rather low luminosity (4 x 10 32 
ergs s _1 ); there is no agreement on whether it belongs to the class of AXPs. There is a claimed 
detection of this source in R and B, but it is unclear whether the observed emission is from RXJ 
0720 or if it is instead due to a background star. If it is the latter, then, for the estimated distance 
of 400 pc, the accretion model of CHN requires the age of this object to be ^ 2 x 10 6 yr for 
the emission to be below the detected values. Since no P has been measured, there is no other 
constraint on the age. If we assume that the optical emission does come from RXJ 0720, then the 
CHN disk model would require an age around 2 x 10 6 yr to produce an emission in roughly the 
same range. The slope of the spectrum, however, is not consistent with the disk model; the model 
predicts that the R band should have higher flux than the B band, while the observations show 
the contrary. Observations in other bands are needed to confirm the spectral sloped]. 



6 Note, incidentally, that also in the case of PRS 0656+14 (Fig. 3), the flux in the B band does not follow the 
trend of the other data points. 
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3.2. Radio pulsars 

In the accretion picture, AXPs and pulsars are drawn from the same population, and thus 
it is reasonable to expect that disks from fallback material should also be found around some 
pulsars. Figure 2 shows the emission spectrum from a disk surrounding a pulsar at various ages. 
The X-ray luminosity of the star at each age has been taken from the cooling curve for a neutron 
star with an iron envelope and B = computed by Heyl & Hernquist (1998) []. In all cases, we 
chose B\2 = 3, and a disk mass equal to M<± = 0.005Mq at an inclination angle of 60°, as for 
AXPs; however, whereas there is a relatively tight constraint on the specific value of the disk mass 
for AXPs, the disk mass could be much lower (or indeed even zero) for radio pulsars. Thus, our 
predictions are less firm. 

We have applied our model to PSR 0656+14. This is an isolated, middle-aged pulsar which 
shows an enhanced emission redward of the V band (Pavlov, Welty & Cordova 1997). The data 
are shown in Figure 3. The soft X-ray part of the spectrum is well fitted by a thermal component, 
originating from the surface of the neutron star. This part of the spectrum looks similar to that 
observed in other pulsars, such as Geminga. As Pavlov et al. notice, however, the excess emission 
in the R and I bands is not very common, and they fitted it with a non-thermal component 
assumed to be produced in the magnetosphere. Here we consider the possibility that, whereas the 
high energy emission comes from the surface of the star (and is well fitted by a Rayleigh-Jeans 
spectrum), the low-energy emission could be produced by a residual disk. Following the notation of 
Pavlov et al. for the Rayleigh-Jeans part, we fit the spectrum as the sum of these two components: 

x 1(T - 4A M , (14) 

where vq is an arbitrary reference frequency (which we took equal to 5.8 x 10 14 Hz), 
g = 3.12 x Kr 31 G ergs cm~ 2 s" 1 Hz -1 , and G = T e (R 10 / d 500 ) 2 . Given the age r = 1.1 x 10 5 yr 
(Pavlov, Stringfellow & Cordova 1996), a distance of 760 pc (Taylor, Manchester & Lyne 1993), an 
extinction of Ay = 0.15 and E(B — V) = 0.05 (Pavlov, Stringfellow &z Cordova 1996), a magnetic 
field B u = 4.7, and an A-ray luminosity of 3.5 x 10 32 (£>/760 P c ) 2 (^/5 X 10" 12 ergs' 1 cm" 2 s" 1 ) 
erg s™ 1 (Cordova et al. 1989, Finley, Ogelman &: Kiziloglu 1992), the only free parameters of 
the disk model are the disk mass and the inclination angle of the disk. Figure 3 shows the 
resulting spectrum corresponding to three choices of these parameters. Notice, however, that the 
observations imply rather low values of M, and the best fit is obtained for an extreme value of the 
inclination angle. 

If the assumption of an accretion disk surrounding this pulsar is correct, then it should 
be possible to observe it directly. Figure 4 shows the predicted spectrum over a wide range of 



7 In the model made by Heyl & Hernquist, various cooling curves, corresponding to different choices of the 
composition of the envelope of the neutron star and strength of its magnetic field, cross around 10 5 yr. Therefore, 
except for either very young or very old systems, our results are not very sensitive to this particular choice. 
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wavelength. The emission in the submillimeter is within the range of sensitivity of SCUBA. 

4. Conclusions 

We have computed the emission from time-dependent accretion disks around neutron stars, 
established by fallback following a supernova explosion. Such disks have been proposed as 
candidates for the energy source of AXPs. We applied our model to predict the emission in the 
optical, infrared and submillimeter bands for the known AXPs. The predicted fluxes are within 
the range of sensitivity of current instruments. We need to emphasize, however, that whereas our 
specific estimates have been made using the disk models of CHN, the arguments should apply in 
general to any accretion model for AXPs regardless of the source of the disk mass. Therefore, 
deep observations at long wavelengths can be used to test all variants of the accretion hypothesis 
for AXPs. 

If accretion disks are indeed responsible for powering the X-ray emission from AXPs, then 
they should also be found around typical radio pulsars. We computed the emission for parameters 
typical of radio pulsars and showed that emission from a disk might explain the excess flux 
redwards of the B band observed in the pulsar PSR 0656+014. Sensitive observations at longer 
wavelengths would severely test this interpretation. 

In addition to perhaps exploring the energy source powering AXPs, accretion disks around 
young compact objects may have other observable consequences. Planets have been found 
around at least one pulsar, PSR B1257+12 (Wolszczane & Frail 1992; Wolszczan 1994). Fallback 
debris would appear to be an ideal source of planet-forming material since it would naturally 
be metal-rich if it is drawn from the inner regions of the progenitor star (e.g. Lin et al. 1991). 
In some cases, fallback may be sufficiently intense to yield a black hole remnant rather than a 
neutron star. Recently, a class of radio-quiet X-ray point sources has been identified, consisting of 
approximately ten members, the majority of which are associated with young supernova remnants 
(e.g. Brazier &: Johnston 1999). Time variability in the emission from the source in SNR RCW 
103, for example, has motivated the possibility that this source may be an accreting black hole 
rather than a cooling neutron star (e.g. Gotthelf, Petre & Vasisht 1999). Qualitatively, accreting 
disks established by fallback around black holes should be similar to those examined here. Thus, 
deep observations at optical, infrared and submillimeter wavelengths may also reveal the true 
nature of radio quiet point sources in supernova remnants. 

This work was supported in part by NSF grant AST 980686. 
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Table 1. — Predicted flux at frequencies v x = 2000 GHz, v 2 = 1.36 x 10 14 Hz (K band), 
u 3 = 5.4 x 10 14 Hz (V band) for some of the known AXPs. The model parameters adopted here are 
the same as in Figure 1. References for these objects can be found in Corbet et al. (1), Sugizaki 
et al. (2), Wilson et al. (3), Osterbroek et al. (4), Gotthelf & Vasisht (5). 



-14- 




Fig. 1. — Emission spectrum from an accretion disk around an AXP at different ages. The disk 
mass is Ml = 0.005M© and the magnetic field strength is B = 8 x 10 12 G, similar to values 
considered in the accretion model of CHN. The irradiating X-ray flux is assumed to derive from 
accretion onto the neutron star. The four curves correspond to ages of r = 10 3 yr (solid line), 
r = 10 4 yr (dotted line), r = 10 5 yr (dashed line), r = 10 6 yr (dotted-dashed line). The inclination 
angle of the disk is assumed to be i = 60°. 
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Fig. 2. — Emission spectrum from an accretion disk around a normal pulsar at different ages. Here 
the disk mass has the same value as in Fig. 1 for AXPs, but the magnetic field strength is smaller, 
B = 3 x 10 12 G, making accretion impossible according to the model of CHN. The irradiating 
X-ray flux is assumed to be generated by cooling of the neutron star. The four curves correspond 
to ages of r = 10 3 yr (solid line), r = 10 4 yr (dotted line), r = 10 5 yr (dashed line), r = 10 6 yr 
(dotted-dashed line). The inclination angle of the disk is assumed to be i = 60°. 
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Fig. 3. — The excess flux redward of the B band observed for the pulsar PSR 0656+14 is consistent 
with emission from an accretion disk. Age, irradiating X-ray flux, dust extinction and distance are 
taken from observations. The various lines correspond to a model with M® = O.OOO15M , i = 60° 
(dashed line), = 0.00023M o , i = 80° (dotted line), = O.OOO4M , i = 87° (solid line). The 
emission at high energies is assumed to come from the surface of the neutron star. 
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Fig. 4. — Emission spectrum from PSR 0656+14 if an accretion disk is responsible for the excess 
emission observed redward of the B band. The various lines correspond to the same models as in 
Figure 3. 



